Cosmological consequences of a possible A-dark matter interaction 
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We propose a general class of interacting models in which the interaction between the CDM com- 
ponent and A is parameterized by an arbitrary function of the cosmic scale factor e(a). Differently 
from other dynamical A scenarios in which the final stage of cosmic expansion is a de Sitter phase, 
we find solutions of transient acceleration, in which the A-dark matter interaction will drive the 
Universe to a new dark matter-dominated era in the future. We investigate some cosmological con- 
sequences of this model and discuss some constraints on its parameters from current SNe la, BAO, 
CMB and Ho data. 
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I. INTRODUCTION 

According to current observational results, a mysteri- 
ous field, named dark energy, which accounts for ~ 70% 
of the cosmic composition, is governing the late time dy- 
namics of the Universe. However, although fundamental 
to our understanding of the Universe, several important 
questions about the nature of this dark energy compo- 
nent and its role in the cosmic dynamics remain unan- 
swered (see, e.g., [l| for some recent reviews). 

Among many possible candidates, perhaps the sim- 
plest explanation for current observations is that the un- 
clumped form of energy density corresponds to a pos- 
itive cosmological constant A, whose presence modifies 
the Einstein field equations to 



(1) 



where G'"' is the Einstein tensor, T'"' is the energy- 
momentum tensor of matter fields and CDM particles, 
and X = SttG is the Einstein's constant (throughout this 
paper we work in units where the speed of light c = 1). 

From the observational point of view, it is well known 
that flat models with a very small cosmological term 
(pa ^ 10"'*'' GeV'') are in good agreement with almost 
all sets of cosmological observations, which makes them 
an excellent description of the observed Universe. From 
the theoretical viewpoint, however, at least two problems 
still remain. First, and possibly the most serious one is 
the unsettled situation in the particle physics/cosmology 
interface (the so-called cosmological constant problem 
(CCP) 3), in which the cosmological upper bound dif- 
fers from theoretical expectations (pA ~ 10^* GeV'') by 
more than 100 orders of magnitude. The second is that, 
although a very small (but non-zero) value for A could 
conceivably be explained by some unknown physical sym- 
metry being broken by a small amount, one should be 
able to explain not only why it is so small but also why 



it is exactly the right value that is just beginning to dom- 
inate the energy density of the Universe now. Since both 
components (dark matter and dark energy) are usually 
assumed to be independent and, therefore, scale in differ- 
ent ways, this would require an unbelievable coincidence, 
the so-called coincidence problem (CP). 

A phenomenological attempt at alleviating the CP 
problem is allowing the dark matter and dark energy to 
interact (in the case of a A-dark matter interaction, A is 
necessarily a time-dependent quantity, which is the basic 
idea behind the decaying A models whose aim is to solve 
or alleviate the CCP)*. Cosmological scenarios with a dy- 
namical A term were independentlyproposed about two 
decades ago in Ref. Q (see also [1, i, S 0, i, i, [3) 
whereas models of couped quintessence (in which the 
dark energy is represented by a scalar field </> or by a 
smooth component parameterized by an equation of state 
Pde = wPjWj with w < Q) have been investigated more 
recently [li [^, la]. In both cases, however, the ab- 
sence of a natural guidance from fundamental physics on 
a possible interacting or coupling term between the two 
dark components leads most of the current investigations 
discussed in the literature to a phenomenological level. 

In this regard, a still phenomenological but very inter- 
esting step toward a more realistic interacting or coupling 
law was recently discussed in Ref. Q, in which the time 
dependence of the A term is deduced from its effect on 
the CDM evolution. Such a coupling is similar to the one 
obtained in Ref. [H from arguments based on renormal- 
ization group and seems to be very general, having many 
of the previous attempts as a particular case. Thermo- 
dynamical considerations for this class of dynamical A 
models [s^] showed that the interacting parameter e must 
be positive, which means that the energy transfer be- 
tween A and the dark matter field is such that the latter 
always gains energy from the former, and not the other 
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^ Strictly speaking, in the context of classical general relativity any 
additional A-type term that varies in space or time should be 
thought of as a new time- varying field and not as a cosmological 
constant. Here, however, we adopt the usual nomenclature of 
time-varying or dynamical A models. 
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way around. In Refs. @, a scalar field description for this 
class of dynamical A models was investigated and found 
to be well represented by a coupled double exponential 
potential of the type V{^) oc exp (A0) + exp {—X(f>). In 
Ref. 12], a coupled quintessence model based on the 
above arguments (for which the equation-of-state param- 
eter w —1) as well as some observational constraints 
on the interacting term were discussed. 

An important aspect worth emphasizing is that in the 
above analyses the interacting parameter e has been con- 
sidered constant over the cosmic evolution whereas in a 
more realistic case it must be a time-dependent quantity. 
Our goal in this paper is therefore to go a step further 
in the above description and extend the arguments of 
Refs. 0,11, Q to a physically more realistic case in which 
e is a function of time. We restrict the present analysis to 
coupled quintessence models in which = — 1, which is 
mathematically equivalent to dynamical A scenarios. In 
terms of the interacting function e(a), we discuss the dy- 
namical behavior of this class of models and find viable 
cosmological solutions for a subset of values of e(a). The 
possibility of transient accelerating solutions in which the 
Universe will experience a future dark matter-dominated 
phase is also explored. We also carry out a joint statis- 
tical analysis with recent observations of SNe la, BAO, 
CMB and Hq to check the observational viability of this 
general class of interacting A-dark matter models. 



II. THE MODEL: BASIC EQUATIONS 

From Eq. ([T]) the Bianchi identities imply that the cou- 
pling between a dynamical A term and CDM particles 
must be of the type^ 



\ X 



or, equivalently, 



Pdm + i — Pdm — — PA , 

a 



(2) 



(3) 



where pdm and pA are the energy densities of CDM and A, 
respectively, and T^^'^ = pdmU^u^ stands for the energy- 
momentum tensor of the CDM field. As usual, the dot 
sign denotes derivative with respect to the time. 

The A-CDM interaction implies that the energy den- 
sity of this latter component must dilute at a different 
rate compared to its standard evolution, pdm oc , 
where a is the cosmological scale factor. Thus, the de- 
viation from the standard dilution may be characterized 



by the function e(a), such that 



dm — PdmfiO, 



~3-|-e(a) 



(4) 



where we have set the present-day value of the cosmo- 
logical scale factor oq = 1. Since the other matter fields 
(radiation and baryons) are separately conserved, Eqs. 
© and dH) provide 

/•I e(5) + ae'ln(a) ^ 
PA = Prfm,o j -4_^(^) aa -I- X , (5) 

where a prime denotes derivative with respect to the scale 
factor and X is an integration constant. 

Neglecting the radiation contribution, the Friedmann 
equation for this dynamical A cosmology can be rewritten 
as 

H = Ha [nbs)a-^ + f^rf,„^o'y5(a) + f^x.o] ^'^ , (6) 

where H(z), Qb.o and ^dm.o are, respectively, the Hub- 
ble parameter, and the baryons and CDM present-day 
density parameters. The function ip{a) is written as 

^ia) = + / <a) + a.-ln(a) ^, 

and r2x,o stands for the present-day relative contribution 
of the constant X to the expansion rate. 

In order to proceed further, we must assume an appro- 
priated relation for e(a). Certainly, among many possible 
functional forms, a very simple choice is 



e{a) 



(8) 



where eg and ^ may, in principle, take negative and pos- 
itive values. With the above expression, Eq. ([5]) can be 
rewritten as 

/■I [1 + ln(5«)] ^ 
PA = Pmoeo J ^^^z^^^da + X. (9) 

Note that, in the absence of a coupling with the CDM 
component, i.e., eo = 0, we may identify X = pAo (the 
current value of the vacuum contribution) and the stan- 
dard ACDM scenario is fully recovered. Note also that, 
for ^ = and eo 7^ 0, the above expressions reduce to the 
dynamical A scenario recently discussed in Refs. 0,[1,[1], 
whose vacuum energy density is given by 



PA-^^^^a-3+^°+X 
3 - eo 



(10) 



In our analysis, we will consider only interaction between A and 
CDM particles. For a discussion about bounds on the interaction 
with conventional matter from local gravity experiments, see Il4( . 
For constraints from primordial nucleosynthesis, see also g. 



III. COSMIC EVOLUTION 

The time evolution of the density parameters r2b(a), 
f^dm(«) and riA(a) (the relative contribution of pA [Eq. 



3 




-8 -6 -4 -2 2 4 -8 -6 -4 -2 2 4 -8 -6 -4 -2 2 

log(a) log(a) log{a) 



FIG. 1: Evolution of the density parameters Qj (j = b,dm,X) as a function of log(a) for some selected combinations of 
eo = ±0.1 and ^ = ±0.2 and A ~ 17.27, B = 0.17 and D ~ 3. 
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FIG. 2: The same as in Fig. 1 for combinations of eo = ±0.1 and f — ±1.0. Note that for eo > and large positive values of ^ 
(Panel Ic), the A-dark matter interaction will drive the Universe to an eternal deceleration instead of the usual de Sitter phase. 



|9] to the expansion rate) can be derived by combining 
Eqs. Q-®. They read 



A + a-3 + B-V(a) ' 

a-3+e(a) 

^ D + Ba-3 + (^(a) ' 

D + tp{a) - a-3+^(") 
D + Ba-3 + ip{a) ' 



(11a) 



(lib) 



(He) 



where A = V,xfi/^b,o, B = flb,o/^dmfi and D = 

Figures 1 and 2 show the evolution of the density pa- 
rameters rij (j = 6, dm, A) with the logarithm of the scale 
factor log(a) [Eqs. (fTTa|) - (|TTc)) ] for values of A ~ 17.27, 
B = 0.17 and D ~ 3 (corresponding to flb^ — 0.0416 and 
^dm,o — 0.24) and eo = ±0.1. Two symmetric values 
of i.e., C = ±0.2 (Fig. 1) and ^ = ±1.0 (Fig. 2) are 
considered. Although currently accelerated (and, there- 
fore, possibly in agreement with SNe la data), models 



with eg > and negative values of ^ (Figs, la and 2a) 
fail to reproduce the past dark matter-dominated epoch, 
whose existence is fundamental for the structure forma- 
tion process to take place. In both cases, the dark energy 
and dark matter densities vanish at high-z and the Uni- 
verse is fully dominated by the baryons. Note that the 
same is not true when both ^ and eo take negative values 
(Fig. Ic). In this case, the negative signs compensate 
each other so that the dark matter dominates the past 
evolution of the Universe while the baryonic and dark 
energy densities vanish. 

Regardless of the sign of eg, well-behaved scenarios are 
obtained when ^ takes positive values (Figs, lb and 2b). 
Note that in these cases a mix of baryons (< 20%) and 
dark matter (> 80%) dominates the past evolution of 
the Universe whereas the dark energy is always the dom- 
inant component from a value of a* < 1 on. A very 
interesting and completely different future cosmic evolu- 
tion is obtained when eo > and the parameter ^ takes 
large positive values (> 0.8). This is shown in Fig. 2c 
for ^ = 1.0 and eo = 0.1. Note that, besides having a 
well-behaved past evolution and being currently accel- 
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log(a) 



FIG. 3: Deceleration parameter as a function of log(a) for 
some selected values of eo and In agreement with the results 
shown in Figs. 1 and 2, note that for large positive values of 
^ the Universe will experience a new matter-dominated era in 
the future, when a —^ oo. 
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FIG. 4: The effective EoS w^ffi^a) = pT{a) / pT^a) as a func- 
tion of log(a) for some selected values of eo and ^. Depending 
on the combination of values for eo and ^ We// rnay bahave 
as freezing or hybrid. 



erating, the cosmic acceleration will eventually stop at 
some value of a >> 1 (when the dark energy becomes 
sub-dominant) and the Universe will experience a new 
matter-dominated era in the future, when a ^ oo. This 
kind of dynamic behavior is not found in most of the dy- 
namical A or coupled quintessence models discussed in 
the literature, being essentially a feature of the so-called 
thawing [l^ and hybrid potentials, which in turn 
seems to be in good agreement with some requirements 
of String or M theories, as discussed in Ref. [17| (see also 

To better visualize this transient acceleration phe- 
nomenon, we derive the deceleration parameter q = 
—ad /a?', given by 



q{a) 



3 n^g-^ + rirf,„,oa^(°)-3 



-1, 



(12) 



and shown in Fig. 3 as a function of log(a) for some 
selected values of ^ and eo- Note that for large posi- 
tive values of ^ the Universe was matter-dominated in 
the past [q{a) -^1/2 for a « 1], switched to a long 
period of cosmic acceleration at aacc < 1 but will even- 
tually decelerate again at some a^ec > 1 (see also |T3| for 
a discussion on quintessence and brane-world models of 
transient acceleration). 



For the sake of completeness, we also show in Fig. 4 the 
effective equation-of-state (EoS) parameter [weff{a) = 
PT{a)/pT{a)] 



Weffia) 



(13) 



as a function of log(a) for some combinations of ^ and eg. 
Note that, while the value of ^ determines the general 
behavior of w{a), the interacting parameter eo is directly 
related to the duration of the accelerating phase. Note 
also that, although presenting many different behaviors, 
clearly a very interesting one is provided by values of 
^ ~ 1, in which w{a) behaves initially as freezing over 
all the past cosmic evolution, is approaching —1 today 
(in agreement with current observations), will become 
thawing in the near future and will behave as such over 
the entire future evolution of the Universe. This freez- 
ing/thawing or hybrid behavior, originally discussed in 
Ref. (l6j . is particularly interesting because, in principle, 
it could reconcile the slight preference of the SNe la and 
large scale structure data for freezing EoS pointed out 
in Refs. 0, [U, (which in turn leads to an eternally 
accelerating Universe) with the String/M-theory require- 
ments discussed in Ref. 1 1711 . 



IV. OBSERVATIONAL ANALYSIS 



^ The argument presented in Ref. [TtII is that an eternally accel- 
erating universe, a rather generic feature of many quintessence 
scenarios (including the standard ACDM model), seems not to 
be in agreement with String/M-theory predictions, since it is 
endowed with a cosmological event horizon which prevents the 
construction of a conventional S-matrix describing particle inter- 
actions. 



From now on we will discuss more quantitatively the 
observational viability of the class of interacting scenarios 
discussed above. To this end we perform a joint analy- 
sis involving current SNe la, BAO, CMB and Hq data. 
Since we are particularly interested in bounds on the pa- 
rameters eo and ^ we fix fif, q = 0.0416 from WMAP 
results (which is also in good agreement with the 
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FIG. 5: The results of our statistical analyses, a) Contours of in the plane eo — These contours are drawn for Ax^ = 2.30 
and 6.17. Note that almost the entire interval of negative values of ^ are ruled out at 2a unless eo — (ACDM model), b) Ax^ 



as a function of the interacting parameter eo. From this analysis, we find e : 



-0.08 



+0.18+0.33 



at 1 and 2cr levels, c) fid 



eo 



space for a constant interacting parameter e = eo. At Icr level eo is restricted to the interval eo = —0.03 ± 0.03. 



bounds on the baryonic component derived from primor- 
dial nucleosynthesis and consider the recent deter- 
mination of the Hubble parameter Hq = 74.2 ± 4.8 [1^ 
in conjunction with the CMB constraint fldm.oh^ = 
0.109 ± 0.006 [23]. 

We use one of the most recent SNe la compila- 
tion, the so-called Union sample compiled in Ref. [2^ 
vifhich includes recent large samples from SNLS and 
ESSENCE [28| surveys, older data sets and the recently 
extended data set of distant supernovae observed with 
the Hubble Space Telescope. The total compilation 
amounts to 414 SNe la events, which was reduced to 307 
data points after selection cuts. 

We also use the distance ratio from zbao = 0.35 to 
zls — 1089, as measured by the Sloan Digital Sky Sur- 
vey (SDSS), Rbao/ls = 0.0979 ± 0.0036 0. Here this 
quantity is given by 



R 



BAO/LS 



-B(^BAO) 



1/3 



r2/3(zBAo) 



rizLs) 



(14) 



where r(z) = dz' / E[z') is the comoving distance. 
Note that, although dark energy does not dominate early 
(see Figs. 1 and 2), we do not use in our analysis either 
the BAO {A) or the CMB shift {TV} parameters since 
these quantities use the approximation that the sound 
travel distance is oc 1/ -^/Om, which is not always true for 
interacting dark matter/energy models due to the process 
of energy transfer (dark matter creation/annihilation [soj 
or varying mass particles (3l| ) between these components 
(see, e.g., [l^ for a discussion). Note also that the above 
quantity seems to be slightly more precise than the A 
parameter since the scatter induced by uncertainties on 
flmh^ cancels out in the ratio. In our analysis, therefore, 
we minimize the function = xInc + Xrbao/ls' ^^i*^^ 
takes into account both the SNe la and BAO / CMB data 



discussed above (we refer the reader to Refs. [32] for more 
on analyses involving different data sets). 

The results of our statistical analyses are displayed in 
Fig. 5. Figure 5a shows confidence contours at 68.3% 
and 95.4% in the parametric space eo ~ C that arise from 
the joint analysis described above. As expected, we note 
that the current observational bounds on ^ are quite weak 
since it appears as a power of the scale factor in the en- 
ergy density [Eqs. ^ and Q]. Note also that very large 
values of the interaction parameter eo are completely ex- 
cluded regardless the value of ^, and that the same is also 
true for almost the entire interval of negative values of ^ 
unless eo ~ which, irrespective of the dimensionless pa- 
rameter f , behaves very similarly to the standard ACDM 
model. 

To better visualize the constraints on eo, in Fig. 5b 
we plot the plane eo — Ax^. From this analysis we 



find eo 



-0.08 



+0.18 
-0.28 



at la level (xLn/^ = 1-02) 



which means that both negative and positive values for 
the interacting parameter are allowed. Physically, this 
amounts to saying that not only is an energy flow from 
dark energy to dark matter (eo > 0) observationally al- 
lowed, but so is a flow from dark matter to dark energy 
(eo < 0) [see Eq. (|3])]. If we take in account the thermo- 
dynamical constraint derived in Ref. Q, i.e., eo > 0, we 
find eo = O.Otol at Icr level. 

For the sake of completeness, we also show in Fig. 5c 
the space fldm,o — eo for the case in which the interacting 
term is constant (i.e., ^ = 0) 0, H, Q. Although phys- 
ically more realistic and producing viable cosmic histo- 
ries as shown in Figs. 1 and 2, clearly the introduction 
of a time-dependence on e(a) (quantified by the param- 
eter ^) weakens the constraining power of the analysis. 
For the case e — const., the interacting parameter eo is 
more tightly bounded, i.e., eo = -0.03 ± 0.06 at 95.4% 
(C.L.). If we consider the constraint eo > Q, we find 
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eo = 0.00l°:°^+;; o6 (68.3% and 95.4% C.L.). 

V. FINAL REMARKS 

We have discussed some cosmological consequences of 
an alternative mechanism of cosmic acceleration based 
on a general class of A-CDM interacting scenarios whose 
interaction term e is deduced from the effect of the dark 
energy on the CDM expansion rate. We have gone a step 
further in the above description and extended the argu- 
ments of Refs. 0, to a more realistic case in which e 
is a function of time [see Eq. ^]. The resulting expres- 
sions for the model are parameterized by the dimension- 
less parameters eo and ^ and have many of the previous 
phenomenological approaches as a particular case. 

We have also investigated the dynamical behavior of 
these scenarios and found a number of viable cosmolog- 
ical solutions for a subset of the parameters eo smd ^ 
(Figs. 1 and 2). In particular, for large positive values 
of ^ (> 0.8) and eo > we have found solutions of tran- 
sient acceleration, in which the A-dark matter interaction 
will drive the Universe to a new matter-dominated era in 
the future. As mentioned earlier, this kind of solution 
seems to be in agreement with theoretical constraints 
from String/M theories on the quintessence potential 
or, equivalently, on the dark energy equation-of- 
state w, as discussed in Ref. [l7j . 

From the observational point of view, we have inves- 



tigated the current bounds on the parameterization ([S]) 
from recent data of SNe la (Union sample) , the distance 
ratio from baryon acoustic oscillation at zbao = 0-35 and 
CMB decoupling at zls = 1089, and Hq estimates. We 
have shown that negative and positive values for the in- 
teracting parameter eo are observationally allowed, which 
means that both an energy flow from dark energy to dark 
matter as well as a flow from dark matter to dark energy 
are possible. For the parameter ^ we have also found 
that positive values are largely favoured over negative 
ones. This includes all the well-behaved cases shown in 
Figs. 1 and 2. 

Finally, it is worth emphasizing that we have restricted 
the present analysis to coupled quintessence models in 
which = — 1 (dynamical A models) whereas a full 
treatment of the dark matter-dark energy interaction 
must also take into account the role of the dark energy 
equation-of-state in the process. Some theoretical and 
observational consequences of a w-CDM interacting sce- 
nario with a time-dependent coupling term, as well as a 
scalar field description for this class of models will appear 
in a forthcoming communication (33| . 
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